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ON  THE  FREE  FLIGHT  MOTION  OF  MISSILES  HAVING  SLIGHT 
CONFIGURATIONAL  ASYMMETRIES 


John  D.  Nicola ides 
Bureau  of  Ordnance 
Department  of  the  Navy 

ABSTRACT 

The  theories  for  the  free  flight  motion  of  missiles,  as 
generally  considered  toy  the  aero-dynamicist  and  the  ballistician, 
are  combined  to  yield  a single  theory  for  the  basically  sym- 
metrical missile.  The  force  and  moment  system  contains  not  only 
the  usual  aerodynamic  forces  and  moments  but  also  the  effecto  of 
slight  configurational  asymmetries  and  the  effects  of  rolling 
velocity. 

The  theory  yields  the  condition  for  the  dynamic  stability 
of  both  statically  stable  and  statically  unstable  missiles,  and 
also  predicts  that  the  presence  of  configurational  asymmetries 
together  with  rolling  velocity  may  result  in  "resonance  insta- 
bility." 

Numerical  integrations  of  the  differential  equations  for 
the  pitching  and  yawing  motion  are  carried  out  for  three  varia- 
tions in  the  rolling  motion.  The  results  indicate  that  the 
rapidity  of  passage  through  the  resonance  region  is  a significant 
factor  affecting  the  magnitude  of  the  pitch  and  yaw  of  the 


missile . 


Two  models  of  a simple  arrow  type  missile  having  control 
surface  deflection  are  gun-launched  at  supersonic  velocity  in 


the  Aberdeen  Spark  Photography  Range  and  the  free  flight 
pitching  and  yaving  motion  and  the  transverse  displacement  are 
measured.  The  tricyclic  theory  is  fitted  to  the  experimental 
data.  The  results  indicate  that  the  theory  accurately  repre- 
sents the  actual  motion  of  the  two  models  and  that  the  as- 
sociated static  and  dynamic  aerodynamic  derivatives  are 
accurately  determined. 
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SYMBOLS 


b fin  or  wing  span 

c fin  or  ving  chord 

Cij  aerodynamic  derivatives  (See  Fig.  2) 

-‘u  -(hr 

d diameter  of  missile  body 

I transverse  moments  of  inertia 


axial  moment  of  inertia 


mass  of  missile 
air  density 

total  velocity  of  missile  in  space 

■XT  J 

reference  area  for  aerodynamic  forces  and  moments , generally  ft  d~ 


■‘vSS 


T total  kinetic  energy  of  the  missile 

XYZ  orthogonal  axis  system  fixed  in  misBile  and  rotating  with  it 
(See  Fig.  2) 

X,  Y,  Z aerodynamic  forces  along  the  X,  Y,  and  Z axes  respectively 
L,  Mj  N aerodynamic  moments  about  the  X,  Y,  and  Z axes  respectively 
(See  Fig.  2) 

u,  v,  w components  of  the  total  linear  velocity  of  the  missile  in  space 
along  the  X,  Y,  and  Z axes  respectively 


•'  v'.". v.’.-.w.- ‘K'Sk ..-v  ..v.\ 


3:J  HjT 


t 


ct 


f 


angle  of  attack  of  missile  in  XYZ  system 

angle  of  sideslip  of  missile  in  XYZ  system 


components  of  the  total  angular  velocity  of  the  missile 
in  space  along  the  X,  Y,  and  Z axes  respectively 
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Sg 
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XYZ 


» z 
✓ 
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U,  V,  V 


control  surface  deflection  angle  (Aileron)  or  asymmetry 
about  X axis 

control  surface  deflection  angle  (Elevator)  or  asymmetry 
along  Z axis 

control  surface  deflection  angle  (Rudder)  or  asymmetry 
along  Y axis 

effective  asymmetry  angle 

orthogonal  axis  system  attached  to  missile  (X  axis  lies 
along  axis  of  mass  symmetry  of  missile,  and  Y axis 
constrained  to  lie  in  xy  plane)  (See  Fig.  l) 

VA  V*  ** 

aerodynamic  forces  along  the  X,  Y,  and  Z axeB  respectively 
aerodynamic  moments  about  the  X,  Y,  and  Z axeB  respectively 
(See  Fig.  1) 

components  of  the  total  linear  velocity  of  the  missile  in 
space  along  the  X,  Y,  and  Z axes  respectively 


'Z  .v  angle  of  attack  of  missile,  XYZ  eyetem 
* IX 

^ v 

a s JL  angle  of  yav  of  missile,  XYZ  system 
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xyz 
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components  of  the  total  angular  velocity  of  the  nlaelle 

V*  L«  ^ 

in  apaeo  along  the  X,  and  Z axe*  respectively 
(See  Fig.  1) 

orthogonal  axis  system  fixed  in  spaoe  (See  Fig.  1) 
forces  along  the  X,  Y,  and  z axes  respectively 
moments  about  the  X,  Y,  and  z axes  respectively 
linear  velocities  along  the  x,  y,  and  i axes  respectively 

W>  ^ 

angular  velocities  along  the  X,  Y,  and  s axes  respectively 
• • 

0 *•  i f , total  angular  velocity  of  the  missile  normal 
to  the  axis  of  mass  symmetry 

complex  angle  of  attack  and  yaw 


oomplex  angular  velocity  normal  to  the  missile 


rate  of  change  of  complex  angle  of  attack  and  yaw 


rate  of  change  of  complex  angular  velocity  normal  to 
the  missile 


INTRODUCTION 


Developments  of  theories  for  the  free  flight  motion  of  missiles 

have,  in  general,  proceeded  along  two  separate  paths.  The  aero- 
1 2 

dynamlcist  ' has  been  primarily  concerned  with  aircraft  which, 
although  lacking  rotational  symmetry  and  essentially  non-rolling, 
are  acted  upon  by  a linear  aerodynamic  aystem  vhlch  includes  forces 
and  moments  due  to  control  surface  deflection  and  lags  in  the  flow. 
The  balllsticlan,  3-11  on  the  other  hand,  has  been  primarily  con- 
cerned vlth  rapidly  rolling  symmetrical  missiles  and  has  included 
in  his  treatment  of  the  motion  the  important  gyroscopic  and  Magnus 
effects. 

In  recent  years,  vlth  the  advent  of  the  guided  missile,  the 
rocket,  supersonic  aircraft,  and  modern  finned  ordnance  veapons, 
the  Interests  of  both  groups  have  merged.  Although  limited 
extensions^2"'1'14  of  the  theories  of  both  groups  have  been  undertaken, 
the  essential  differences  remain. 

It  is  one  of  the  purposes  of  this  paper  to  unify  these  ap- 
proaches for  the  class  of  missiles  vhich  are  basically  symmetrical 
and  are  only  slightly  asymmetrical  due  to  control  surface  deflec- 
tion, ving  and/or  tail  incidence,  bent  body,  damaged  or  malaligned 
fin,  etc. 


li 


The  neceeaity  for  considering  thia  union  arlaea  from  the 
failure  of  tha  exlating  theory1"^  to  repreaent  tha  free  flight 
notion  of  winged  and/or  finned  nlaallea  and  ita  failure  to  account 
for  various  phenomena  which  have  haen  generally  experienced  on 
statically  stable  missiles.  Some  of  these  phenomena  are 

(1)  that  non- rolling  statically  stable  missiles  generally 
have  large  dispersion  and  that  rolling  the  missile  re- 
duces the  dlsperaion. 

(2)  that  even  for  generally  veil  performing  missiles  a fev 
go  berserk1'*  yielding  extremely  poor  dispersion  and 
sometimes  tumbling,  and 

(3)  that  peculiarities  In  the  free  flight  motion  seem  to 
occur  when  the  rolling  velocity  and  the  pitching  velocity 
approach  coincidence. 

The  general  procedure  to  be  followed  here  will  be  to  develop 
the  differential  equations  of  motion  for  a missile  having  basic 
trigonal  or  greater  configurational  symmetry1^  and  slight  con- 
figurational asymmetry.  The  aerodynamic  system  will  include  the 
forces  and  moments  generally  considered  by  both  the  aerodynamic let 
and  the  ballistician.  The  important  gyroscopic  terms  resulting 
from  rolling  velocity  are  also  included.  For  simplicity  and 
clarity  the  differential  equations  will  be  solved  for  the  case  of 
constant  axial  velocity  and  constant  rolling  velocity  of  the 
missile.  The  resultant  tricyclic  theory  for  the  free  flight 


pitching  and  yawing  notion  of  the  mist lie  and  the  theory  for  the 
displacement  of  the  center  of  gravity  of  the  missile  will  be 
applied  to  the  experimental  data  obtained  from  the  free  flight 
motion  of  gun-launched  models  tested  in  the  Aberdeen  Spark 
Photography  Range.  The  static  and  dynamic  aerodynamic  de- 

rivatives associated  will  be  obtained.  Finally,  the  differential 
equations  of  motion  will  be  numerically  integrated  for  the  case 
of  varying  rolling  velocity  in  order  to  indicate  the  motione 
obtained. 


THEORY 


The  difficulties  generally  encountered  in  the  formulation  of  a 

theory  of  free-flight  motion  may  be  separated  into  two  groups, 

dynamical  and  aerodynamical.  Herein  the  dynamical  problem  will  be 

approached  by  employing  a modified  Eulerian  axis  system  and  by 

uslng  the  Lagrange  equation  for  formulating  tha  basic  differential 

equations  of  motion.  Mathematical  simplifications  are  Introduced 

by  limiting  the  angular  displacement#  and  angular  velocities 

(except  rolling)  to  small  size.  This  is  the  familiar  dynamical 

S:  approach  to  the  linear*  :•  -vl  motion  of  a "spinning  top"  or  a 

£l  42 

"gyroscopic  pendulum"  ' . The  aerodynamical  problem  will  be 

approached  in  the  manner  and  nomenclature  of  the  aerodynamic 1st. 

A linear  force  and  moment  nyetem  is  assumed  and  the  symmetry 
arguments  of  the  balllstlcian  employed.  In  general,  then,  these 
two  basic  problems  will  be  approached  separately.  Their  resolu- 
tions will  then  be  combined  and  the  fundamental  differential 
equations  of  motion  will  be  obtained. 

The  differential  equations  of  motion  will  be  solved  for  the 
case  of  constant  flight  velocity  and  constant  rolling  velocity, 
and  the  resulting  expressions  for  the  pitching  and  yawing  motion 
and  for  the  trunsverse  displacement  of  a missile  will  be  discussed 
in  detail.  In  a later  section  the  differential  equations  of  motion 
will  be  investigated  numerically  for  the  case  of  constant  flight 
velocity  but  varying  rolling  velocity. 


.4 


DYNAMICAL  SYSTEM 


The  coordinate  systems  illustrated  in  Fig.  1 vill  he  ueed  in 
considering  the  free  flight  motion  of  a missile.  The  xya  system 
is  orthogonal  and  fixed  in  space.  The  XYZ  system  is  orthogonal 
and  is  pitching  and  yaving  hut  not  rolling  vith  the  missile. 

(The  X axis  lies  along  the  axis  of  mass  symmetry  of  the  missile  and 
the  Y axis  is  constrained  to  lie  in  the  xy  plane).  The  angular 
velocity  of  the  missile  is  given  hy  the  components  p,  q,  and  r in 

V*  * w* 

the  XYZ  system* 

The  coordinates  of  the  dynamical  system  are  taken  as 

(1)  X*  y,  z,  tha  components  of  the  linear  velocity  of  the 
center  of  gravity  of  the  missile;  and 

• B » 

(2)  'p  , the  components  of  the  angular  velocity  of 


the  ml salle  shout  lta  center  of  gravity,  (it  should 
he  noted  that  these  components  are  in  a moving  non- 
orthogonal  modified  Eulerlan  axis  system. ) 


The  total  kinetic  energy  of  the  missile  in  free  flight  is  thus 


given  by 


1 Ti.'2*  1 2 _ 1 T,S  2 . 1 _ i 2 1 _v2*  1 

* *7j“  P * * "5"  Izr  * ® x * ® y * 


It  la  asaumed  that  the  missile  has  symmetrical  mass  distribution, 


IS  ■ I 


and  from  Fig.  1,  it  Is  seen  that 
P s i - £ tin  9 

4 s 0 

r s f cos  9 


Thus  the  total  kinetic  energy  of  the  missile  may  be  written  as 


T ■ Ix  jf2  - Ix  0 ^ sin  9 * | Ix  [j?  sin2  9 • | I & 2 ♦ 
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2 . 1 


axe 


— • 2 i • p 

* 5°5 


Nov  considering  the  Lagrange  equation2 3 for  this  dynamical  system 


d / c»  T \ - t - 
« Sir J 3 4r 

where  qr  z coordinates  of  the  dynamical  system 

(yi,  9,  f , x,  y,  z) 


Qj.  s Generalized  Force  (i.e.,  force  or  moment  tending  to 
change  the  particular  coordinate). 


M*-* 


Performing  the  indicated  ope rat loot,  assuming  that  0,  ^ , ©>  p 

are  email  quantities  and  that  their  products  may  be  neglected,  yields 

• • • “ 

i jfZ  - i ix  0 -C2  - ^ xl 

(9) 

(10) 

(11) 

I l 

^*2  - q * ifJ  , total  angular  velocity  of  the  miesile  normal  to 

the  axle  of  maee  lymoetry  (12) 

S ; y » 1 i | traneverse  displacement  of  the  missile  (13) 

Eqs.  (8)  - (ll)  are  then  the  basic  differential  equations  of  motion  of  the 
missile.  They  may  be  completed  once  the  Generalized  Forces  tending  to  change 
the  coordinates  of  the  dynamical  system  are  knovn.  These  Generalized  Forces 
are  derived  in  the  following  section  from  the  aerodynamic  forces  and 
momenta  which  act  on  a missile  in  free  flight. 


^ (IX  0 - IX  f .in  «)  . 40 

eS 

m S ■ Qg 

■ x ■ St 

where  complex  variables  have  been  Introduced  in  defining 
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k 
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The  coordinate  ayrtea  used  In  considering  the  aerodynaaie  forces  and 
sonants  which  not  on  a nlsslle  in  frt*  flight  will  be  tha  standard  I.A.C.A. 
aystsa  which  la  orthogonal  and  fizad  to  tha  aiaalla  (i.e.,  rolling  with  tha 
aiaalla,  principle  axla).  This  ays  tea  ia  dealgnatad  by  ZZZ  and  tha  coapooente 
of  the  linear  and  angular  velocity  are  given  by  u,  v,  w and  p,  q,  r respective- 
ly  (See  fig.  2). 

Aerodynamic  force 

The  total  aerodynaaie  force  which  acta  on  a niaaile  in  free  flight  la 
aeaiaed  to  depend  on  the  linear  velocity  of  tha  alsslle,  the  angular  velocity 
of  the  niaaile,  accaleratlona,  the  denalty  of  tha  air,  the  velocity  of  eound, 
the  alte  and  ahape  of  tha  aiaalla,  and,  for  tha  particular  eaae  under  ecocide  re 
tlon  here,  alight  aayaaMtry  of  tha  configuration  (l.a.,  control  aurfaoe  de- 
flection, wing  and/or  tall  lneldenoe,  bent  fina,  bant  body,  ate.). 

Aasuaing  that  tha  notion  la  confined  to  inflnltaalaala  and  that  tha 
dependence  la  linear  yields 

r.TpP  »*rr.l£(3  .1,  r 

* V »*■  * XK  N ♦ b*  **  * ” • h.  s « 

z S Z^  * ♦ Zq  q * Z^  t Z^  q 

♦ * v pr  * Sp  p P * V P*  • \ $ e 

where  I.  and  Z.  are  the  Stability  Derivatives. 
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LINEAR  VEL.  * 


u 

v 

w 


* V 


FORCE  * 


X*Cx£fV2S 
Y-Cyi  /*V2$ 
Z-CzJ,V2S 


ANGULAR  vel.= 


p 

q 

r 


MOMENT  « 


L*C^  V2  S b 

M'Lli  fVzSc 
N«C„^V2Sb 


FIG.  2 - AERODYNAMIC  SYSTEM 


Multiplying  Eq.  (15)  by  i and  adding  to  Bq.  (lk),  and  replacing  the 
etability  derivatives  by  the  standard  nomenclature  for  the  aerodynamic 
derivatives  yields 


Here  ve  Impose  the  fundamental  assumption  that  the  contribution  of 
the  configurational  asymmetries  is  "Blight"  and  thus  that  not  only  may  it 
be  added  linearly  but  that  the  basic  missile  (i.e.,  the  configuration  of  the 
missile  when  ^ r » <S  j.  ■ 0)  has  trigonal  or  greater  rotation  symmetry  and 

mirror  symmetry.10!  11 » ^ 


Because  of  these  equalities  of  the  aerodynamic  derivatives  for  the 


class  of  configurations  having  basic  symmetry,  an  Aeroballistlc  nomenclature 
is  convenient  and  is  Introduced  as  indicated  above, 


2k 


a** 

&/  • f * 1 °L 


pq*ir  * ft,  * i«6 

s |3  ♦ isc 


* 

X 


• s 


S q * 1 r 


(18) 
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- q t i r 


a 

X 


: q ♦ 1 r 


Substituting  Eq.  (IT)  - (18)  into  Eq.  (1 6)  yields 
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Npot 


(zv)  ^ (iv) 


♦ 1 c* 


> • 1 °»pi  ($)J(W  *[c«Pi  (•)  * 1 * <\iC 


(19) 


vhert 


(s)  C-  5 s 


1 


°*j,s  B * 1 \ Se 


and  CM  Is  real. 


Since  all  the  terms,  except  the  last,  are  independent  of  the  roll 
orientation  of  the  coordinate  system  and  since  the  last  term  represents  an 
asymmetry  which  la  fixed  in  the  missile  and  thus  rolling  wit'  the  missile, 
t\i&  total  normal  force  may  be  expressed  In  terms  of  the  original  dynamical 
XYZ>  system  which  was  not  rolling  with  the  mlaalle  ee 


- .n 

t i z 

♦ b 

C* 

f 4 

c 6l  ♦ d 

IP 

L “v 

' 1 ( CT ) 

*«  1 p 
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Ls. 

($)  * 1 c 

S.  efi.1 
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J Pit 
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Nt 
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c ■JpVZsjOB^  ♦ IC^  ($)J  (fy) 

4 = lpv2B[cV  <It>*1c»4](57>2 

Th«  Generalized  Force,  Qg,  vhlob  le  tending  to  move  the  mleelle  In  the 
yz  plane  la  the  emanation  of  all  the  foroee  acting  in  thie  plane.  This 
includes  not  only  the  aerodynamic  force*  but  also  the  gravitational  force, 
thrust  malalignment  forces,  and  any  others  that  might  be  acting.  However, 
in  order  to  emphasize  the  effects  of  oonfiguratlonal  asymmetries  and  to 
keep  the  treatment  as  elementary  as  possible,  the  Generalized  Force  will 
be  assumed  to  Include  only  the  aerodynamic  foroes.  (it  should  be  noted, 
however,  that  the  addition  of  the  gravitational  force,10  the  thrust  mal- 
alignment forces  and  mass  malalignment  forces2^  present  no  fundamental 
difficulties. ) 

V*  1*  i*  V* 

Since  Y f 1 Z lies  in  the  Y Z plane  which  has  been  assumed  to  make  a 
small  angle  with  the  yz  plane  (l.e.,  Jl  assumed  small),  then  Y ♦ i Z may 
be  taken  as  equal  to  V * ■ * • 1 * 

Aerodynamic  Moment 

The  specification  of  the  aerodynamic  moment  acting  normal  to  the 
missile  in  free  flight  proceeds  in  the  same  manner  as  the  specification 
of  the  normal  aerodynamic  force.  Accordingly  M and  N are  assumed  to  be 
linear  function!  of  oC  , A , q,  r,  <*  , a , q,  r,  <$,.  and  as 


M . ♦ M.  q ♦ M . o(  * ML  q 
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»{S  f * % r ♦ N(i  (3  ♦ H*  r 


N„v  P*  ♦ N p*  f N . p*  ♦ N 


Multiplying  Bq.  (24)  by  1 and  adding  to  gq.  (23),  and  replacing  the 
stability  derivatives  by  the  aerodynamic  derivatives  yields 


However,  from  symmetry  considerations  it  follows  that 


Thus  the  aerodynamic  moment  "becomes 
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This  aerodynamic  moment  written  in  terms  of  the  non-rolling  system 


le  given  by 


I*  — » — • T“  T*  r i \ pot 

H(ll  s*^  ♦BJ  »C^  ^ 4 J ( 

Since  the  moment  normal  to  the  missile's  c,ie,  is  Assumed  to  con- 

✓ » ***X 

tain  only  the  aerodynamic  moments  and  e inee  K «■  i N is  equal  to  U n we  have 


Q-  * K.  •>  i K 
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DIFFERENTIAL  EQUATIONS  OF  MOTION 


General  Differential  Equations  of  Motion 


With  the  specifications  of  the  aerodynamic  forces  and  momenta  acting  on 
the  missile  in  free  flight  in  terms  of  the  Generalized  Forces,  the  differ- 
ential equations  of  motion  may  he  completed  aB 

♦ • *-  - m r 

l ^ 1 £1  • kC*C  • % ♦ C * D ^ ♦ E e1-^  pdt  (3 

j (Ix  i - Ix  ein  ®)  ■ y (2 

wi  W4  **^  f 

WS  • taC  * b^  ♦ e o&  ♦ d^  *JL  &t,  e1  ' ^ (: 

W#  « Q*  (: 

and  under  the  basic  assumption  of  small  angles  and  slowly  changing  angles 

-£■  • 

it  is  seen  from  Sqs.  (4),  (5),  (12),  and  (10)  that  ^ - jfl.  • (: 

Differential  Equations  of  Motion  For  Constant  Axial  Velocit 


When  ^ is  zero,  it  follows  from  Eq.  (35)  and  the  assumption  of 
small  angles  that  the  axial  velocity  of  the  missile  is  a constant.  The 
sum  of  the  three  small  angles,  namely,  (a)  the  angle  between  the  missile's 
axis  and  the  total  velocity,  (b)  the  angle  between  the  total  velocity  and 
the  fixed  axis,  and  (c)  the  angle  between  the  fixed  x axis  and  the  missile's 
axis,  is  zero. 

Differentiating  this  sum  yields 

0 — ♦ 

S : V ( - i XI  * Oi  ) (3 

Neglecting  products  of  the  aerodynamic  derivatives  in  comparison  with 
the  derivatives  themselves,  Eqs.  (32),  (34)  and  (37)  may  bo  combined  to 


■.vwv/Cv  v, : 
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yield  the  differential  equation  for  the  pitching  and  yawing  motion. 


(33)  and  (30)  will  then  he  considered  in  the  following  sectior. 
First,  the  solution  will  he  given  for  the  oase  of  constant  rolling 
velocity  (i.e.,  ■ 0)  and,  then,  in  a later  section  the  case  of 

varying  rolling  velocity  will  he  numerically  investigated. 


so 


SOLUTION  OF  PITCHING  AND  TAWING  MOTION  FOR  CONSTANT  ROLLING  VELOCITY 

When  Q0  la  zero,  It  follows  from  Eq.  (33)  and  the  assumption  of  small 
angles  that  the  rolling  velocity  of  the  missile  is  a constant.  For  this  case 
the  general  solution  of  the  differential  equation  for  the  pitching  and 
yawing  motion,  Eq.  (38),  is  given  hy 


0</  s K;  e ^ 1 t ♦ Kg  e ^ t ♦ e ^3  t 
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where 
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A physical  representation  of  this  solution  Is  given  by  noting  that  the  motion 
Is  "trieylic";  that  Is  to  say,  the  free  flight  pitching  and  yawing  motion  of 
missiles ' having  slight  configurational  asymmetry  may  he  represented  by  the 
motion 'traced  out  by  three  rotating  vectors.  (See  Tig.  3)*  Rewriting  Eq. 
(40)  as 


^ (Xx  •!«!)*  (>2 

i K.  e ♦ Kg  e 


♦ K 
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ipt 


(46) 


where  Xj.  and  \ 2 are  the  real  parts  and  6)^  and  & 2 ara  the  imaginary  parts 
of  ^i2.  It  is  noted  that  the  real  parts,  \ 1 and  \ 2,  c*Uia  the  and 
Kg  vectors  to  damp  or  expand  and  that  the  Imaginary  parts,  ^ and 
cause  the  vectors  to  rotate.  Since  ip  is  a pure  Imaginary,  the  K3  arm  . 
does  not  change  in  size  but  rotates  at  a constant  angular  velocity  equal 
to  the  steady  rolling  velocity  of  the  missile. 


Before  considering  the  effects  of  configurational  asynmetries  It  is 
helpful  to  review  the  free  flight  pitching  and  yawing  motions  obtained 
for  the  case  of  no  configurational  asymmetries. 
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When  the  asymmetry  la  set  equal  to  zero,  <5C  ■ 0,  the  aolution,  Eq 

ft 

a 

(40),  reduces  to  the  eplcylic  form.3* 
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This  aolution,  like  the  original,  Eq.  (40),  applies  to  both  atatloally 

stable  and  atatloally  unstable  missiles  (l.e.,  to  uiaailea  whose  center 

of  pressure  of  the  normal  force  due  to  angle  of  attack  and  yaw  la  aft  of 

center  of  gravity  of  the  missile  ( -Cm  ) and  to  missiles  whose  center  of 

<k. 

pressure  is  forward  of  the  center  of  gravity  ( ♦ CM  ) . 


Non-Rolling  Missiles.  For  the  case  of  no  rolling  velocity  (p  ■ o) 
the  constant#  are  given  by 

I.  Statically  8table  Missiles  (-  ) 


II.  Statically  Unstable  Missiles  ( t Cu  ) 


For  the  statically  stable  missile  the  vectors  and  Kg,  rotate  in 
opposite  directions  vith  equal  velocity,  and  thus  the  pitching  and  yavlng 
motion  la  given  by  lines,  ellipses,  or  circles  (See  Fig.  4). 


The  general  condition  for  dynamlo  stability  (i.e.#  the  requirement  that 
the  notion  repeat  lteelf  or  damp  out}  ie  that  ^ ^ 2 s 0*  For  the  ease  of  the 

etatically  stable  non- rolling  missile  this  dynamic  stability  condition 
reduces  to 

l\  **'*  \j'lk's  \ ! (,k) 
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since  and  are  negative  and  la  positive. 


For  the  statically  unstable  missile  the  and  Kg  vectors  do  not 
rotate  and  therefore  the  motion  la  a line.  The  motion  is,  however, 


generally  dynamically  unstable  since 


and  thus  one  of  the  arms  vlll  damp  and  the  other  expand. 


Rolling  Missiles.  The  motion  and  the  effects  of  the  aerodynamic 
derivatives  on  the  stability  of  the  rolling  missile  are  more  readily 
discussed  if  the  radical  in  Sqa.  (44)  is  approximated  by  a binomial  expan- 
sion^ and  if 


If2 


Accordingly,  the  real  and  imaginary  parts  of 
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are  given  by 
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subject,  however,  to  the  condition  that 
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For  the  statically  stable  missile  It  Is  seen  from  £q.  (5)  that  the  two 
vectors  rotate  In  the  opposite  direction,  as  In  the  non-rolling  case,  but 
that  the  difference  in  their  rotational  velocities  depends  on  the  rolling 
velocity  of  the  missile.  For  small  rolling  velocity  of  the  missile,  the 
pitching  and  yawing  motion  is  given  by  slowly  rotating  ellipses.  As  the 
rolling  velocity  increases  the  motion  becomes  flower- like  and  finally  for 
large  rolling  velocity  we  have  motions  characteristic  of  the  fast  spinning 
gyroscopic  pendulum.  (Bee  Fig.  5) 


I 

as 
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For  the  statically  unstable  mieeile  it  ie  eecn  from  Eq.  (57)  that  the 
vectors  rotate  la  the  same  direction  (which  ie  aleo  the  eame  direction  ae  the 
rolling)  end  tfcue  give  the  familiar  motione  characteristic  of  a "top”. 

(See  Fig.  6). 

The  effects  of  the  various  aerodynamic  derivatives  on  the  dynamic 
stability  of  a missile  are  more  readily  seen  by  writing  Eq.  (56)  in  the  form* 


(59) 
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2 (stability  factor  of  the  bal 

(60)  (61) 


liatician)^. 

flanges  of  values  for  the  stability  factor  for  the  various  types  of 
missiles  are  given  in  Table  I. 


Statically  Stable  MiBslleB 

S 

< 0 

^<1 

Statically  Unstable  Missiles 

s 

> 0 

t > 1 

Non-Rolling  Missiles 

B 

= 0 

f*  0 

Neutrally  Stable  Rolling  Missiles 

S 

: re 

r*  1 

TABLE  I 

Values  of  B in  the  range  0460  may  not  be  considered  when  using 
Eq.  (59)  due  to  the  condition  imposed  by  Eq.  (58). 


® Eq.  W Tn“similar  form  is  given  by  s f \ t 
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For  statically  stable  rolling  missiles  it  la  seen  front  Bq.  ($9)  that 

'o  and  (L  tend  to  make  the  motion  dynamically  stable,  vhereae  "&!,  and 
"q  •< 

CL  (the  Magnua  moment  coeffiflent  la  generally  taken  aa  positive  for 

V 


the  statically  stable  missile)  have  the  opposite  effect. 


V 


tends  to 


undamp  the  fast  rotating  vector  (nutation)  and  damp  the  slowsr  rotating 
veotor  (precession).  Since  the  size  of  the  term  slowly  increases 

with  rolling  velocity,  there  might  be  limiting  rolling  velocity  beyond  which 
the  missile  oannot  be  dynamically  stable.^ 

For  the  statloally  unstable  missile  the  effects  of  the  aerodynamic 
derivatives  on  the  dynamic  stability  are  not  as  simple.  Prom  Bq.  ($9)  It 
Is  seen  that  C„  tends  to  undamp  the  nutation  and  damp  the  preoeselon, 

y 

whereas  CL  , A , and  CM  (here  the  sign  of  (X.  Is  taken  as  negative, 

Mq  ™ '’pa*-  ”p* 

however,  positive  values  are  not  uncommon)  tend  to  damp  tho  nutation  and 

undamp  the  precession. 


ftynamic  Stability  Criterion 


The  criterion  for  dynamic  stability  is  that 


or  that 


Extracting  the  real  part  of  the  radical,  this  criterion  may  he  replaced 
by  tvo  conditions: 


(1)  when10'26'27 


(2)  vhen  27 


(67) 
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For  the  case  when 


no  dynamic  stability  Is  possible. 

It  should  be  noted  that  although  for  this  case  of  constant  flight 
velocity  and  constant  rolling  velocity  these  conditions  must  be  satisfied, 
for  the  case  of  varying  rolling  velocity  ana  flight  velocity  dynamic  in- 
stability may  be  and  is  tolerated  in  aome  designs  for  short  durations. 

Tricyclic  Motion 

The  presence  of  slight  configurational  asymmetries  has  been  shown  to 
produce  tricyclic  motion  which  differs  from  the  epicyclic  motion  only  in 
the  addition  of  a third  term  and  a modification  in  the  initial  size  of  the 
nutational  and  precesaional  vectors.  Therefore,  the  remarks  of  the  previous 
paragraph  on  the  dynamic  stability  and  the  contribution  of  the  aerodynamic 
coefficients  to  the  motion  still  apply.  It  remains,  however,  to  consider 
the  effect  of  this  additional  third  vector. 

For  the  dynamically  stable  missile  after  the  nutational  and  precesaional 
arms  have  damped  the  third  vector  is  seen  to  represent  the  steady  state 
pitching  and  yawing  of  the  missile. 


Non-Rolling  Missiles.  For  the  statically  and  dynamically  stable  non- 
rolling missile,  the  size  of  the  third  vector  Is  given  by 


vhlch  is  the  familiar  expression  for  the  steady  state  "trim"  of  an  air- 
craft  due  to  elevator  deflection.  * This  non-rolling  trim  angle  and  the 
angle  of  effective  control  surface  deflection  lie  in  the  same  plane  (i.e., 
and  are  parallel).  The  transient  pitching  and  yawing  motion  is  there- 
fore given  by  ellipses,  lines,  or  circles  vhose  centers  are  displaced  from 
the  origin  by  this  angle  of  trim.  (See  Fig.  7)  The  case  of  the  statically 
unstable  non-rolling  miBSile  is  handled  by  the  theory  but  la  hardly  worth 
discussion  since  the  motion  is  generally  dynamically  unstable.  (See  Eq. 
(52).) 

Rolling  Missiles.  The  addition  of  rolling  velocity  produces  profound 
changes  In  the  nature  and  size  of  the  pitching  and  yawing  motion.  The  size 
of  the  third  vector  aa  effected  by  the  rolling  velocity  is  best  Been  by 
conaidering  the  ratio  of  the  rolling  trim  to  the  non- rolling  trim. 
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(71) 


This  ratio  is  called  the  "magnification  factor"  and  typical  value*  aa  a 
function  of  rolling  velocity  are  given  in  Fig*  8 which  were  calculated  for 
the  statically  and  dynamically  stable  missile  used  in  the  experimental  tests 
to  be  reported  in  a later  section.  From  Eq.  (71)  and  from  Fig.  6 it  1b 
apparent  that  a resonance  phenomenon  appears  possible  if  the  rolling  velocity 
approaches  the  rotational  velocity  of  either  the  nutation  vector  or  the 
precession  vector.  However , for  statically  stable  missiles  it  has  been  Bhown 
that  the  preoeesional  vector  and  the  nutational  vector  rotate  in  different 
directions  and  thus  only  an  equality  of  the  angular  velocity  of  the  nutation 
vector  and  the  rolling  velocity  may  be  considered.  Since  for  the  statically 
unstable  missile  both  the  nutational  vector  and  the  precessional  vector  rotate 


FIG  8 - MAGNIFICATION  FACTOR  FOR  TEST 


result  of  these  conditions  It  Is  seen  that,  In  general,  resonance  can  only 
occur  for  the  case  of  the  statically  stable  missile  and,  then,  only  when 
the  angular  velocity  of  the  nutational  vector  equals  the  rolling  velocity 
of  the  missile. 

The  size  of  the  rolling  trim  at  resonance  depends  on  the  rate  of  damping 
(^1,2  If  the  toroplog  10  zero*  the  size  of  the  rolling  trim  ia  infinite. 
However,  since  most  missiles  are  dynamically  stable  the  size  of  the  rolling 
trim  is  finite.  For  missiles  with  good  dynamic  stability  magnification 
factors  of  10  are  not  unusual.  For  missiles  of  marginal  dynamic  stability, 
values  from  50  to  100  are  quite  possible. 

The  effect  of  this  resonance  phenomenon  is  to  produce  large  angles  of 
pitch  and  yaw  with  the  result  that  the  basic  assumptions  of  the  theory  may 
no  longer  be  valid  and  instability  may  result  from  causes  not  considered. 

For  thla  reason  stability  considerations  must  include  "resonanoe  instability"' 
as  well  as  statia  and  dynamlo  stability. 

Orientation  of  Trim  and  Asymmetry.  The  addition  of  rolling  velocity  not 
only  affects  the  size  of  the  rolling  trim  but  also  the  angular  orientation 
between  the  plnne  c'v’*"',''ing  the  effective  control  surface  deflection  and  the 
plane  containing  the  trim  (i.e.,  the  angle  between  K3  and  ).  This  angle 
1b  given  by 


and  typical  values  are  given  in  Fig.  9 for  the  particular  missile  previously 
mentioned. 


SOLUTION  FOR  TBS  TRANSVERSE  DISPLACEMENT  FOR  CONSTANT  ROLLING  VELOCITY 


The  motioh  of  the  center  of  gravity  of  the  missile  may  now  be  considered. 
The  displacement  of  the  missile  in  the  x direction  hat  already  been  assumed 
to  be  one  of  constant  velocity.  Thus,  the  problem  reduces  to  solving  the 
differential  equation  of  motion  for  the  transverse  displacement  of  the  missile 
which  Is  given  by  Eq.  (3*0*  The  function,^  (t)  is  known  from  Eq.  (40 ) , and 
the  function,  Q (t),  may  be  obtained  from  Eq.  (36)  and  Eq.  (37) 


- i <k  (*)  (7b) 

vhers  oC  (t)  is  obtained  by  differentiating  Eq.  (40).  Substituting  these 
functions  into  Eq.  (34)  and  solving  gives  the  solution  for  the  transverse 
displacement  of  the  missile, 
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vhere  and  kcj  are  functions  of  the  initial  conditions  of  flight. 

This  solution  for  the  transverse  displacement  is  seen  to  be  a tricyclic 
motion  plus  a linear  term  and  a constant.  The  constant  term  depends  upon 
the  selection  of  the  coordinate  system.  For  the  case  where  the  origin  is 
the  muzzle  of  the  gun  or  launcher  and  the  x sxiB  lies  along  the  center  line 


of  the  gun  or  launcher,  the  conatant  of  the  linear  term,  le  called  the 
"jump  angle"  and  the  term  Itself  represents  the  line  shout  which  the 
occillatory  motion  (i.e.,  swerving  motion)  takes  place. 

For  the  case  of  the  non-rolling  missile  the  solution  for  the  transverse 
displacement  reduces  to 


In  the  Assign  and  construction  of  most  fin-stabilized  miss  lias,  particu- 
larly ordnance  veapona  (i.e.,  boats,  mortars,  bazooka*,  rockets,  finned 
artillery  shells,  etc.),  configurational  asymastrlaa  are  present  due  to 
various  causes,  some  of  which  are  manufacturing  Inaccuracies,  damage  in 
handling,  damage  in  launching  and  lntantlonal  design*.  These  asymmetries, 
besides  producing  angle  of  incidence  (&%  ,Sg  ),  have  an  equal  probability 
of  producing  differential  angle  of  incidence  ( ).  This  differential  angle 

of  Incidence  oauses  a varying  rolling  velooity  which  for  most  designs  starts 
at  taro  and  approaches  a steady  state  value  vhlch  is  determined  by  an  equality 
of  the  roll  moment  due  to  differential  angle  of  incidence  and  the  roll  moment 
due  to  rolling  velocity.1^'**0 

Si  • a * s * «*> 

The  differential  equation  of  pitching  and  yawing  motion,  Bq.  (36)  has 
been'  integrated  numerically  for  the  following  rolling  motions : 

(1)  aero  rolling  velooity  to  steady  state  rolling  velooity  equal 
to  the  nutation  rate 

(2)  zero  rolling  velocity  to  steady  state  rolling  velocity  equal 
to  1.3  nutation  rata 

(3)  zero  rolling  velocity  to  steady  state  rolling  velocity  equal 
to  3 nutation  rate 

The  resulting  pitching  tnd  yaving  motions  are  plotted  in  Pigs.  10,  11, 
and  12  for  the  missile  used  in  the  experimental  tests,  it  is  sesn  from  the 
figures  that  the  speed  of  passage  through  the  resonance  region  has  a profound 
influence  on  the  magnitude  of  the  pitching  and  yawing  motion. 

< Canting  the  fins,  for  example,  ia  a standard  practice  in  many  ordnance  dedans. 
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COMMENT  ON  THEORY 


In  awaaary,  tha  aulution*  for  tha  fret  flight  motion*  of  a mlaalla 
having  alight  oonfigurational  aaymmatrlaa  and  flying  at  eonatant  veloolty, 
axial  and  rolling,  indicate*  thatt 

(1)  The  pitching  and  yawing  motion  la  Tricyclic. 

(8)  The  tranaveree  diaplacament  la  Tricyelie,  plua  a linear  term 
and  a eonatant*  (For  the  oaae  of  aero  rolling  velocity  the 
third  am  la  replaoed  by  a quadratic  tana*) 

(3)  The  preoeaaion  and  nutation  vector*  rotate  and  changa  alae  aa 
in  the  Bpiaycllo  oaae* 

(4)  The  rotation  and  ohange  in  aiae  of  tha  naw  third  veotor, 
trim,  arlalng  from  the  aaymmetry  depanda  primarily  on  the 
rolling  veloolty  of  the  mlaalla. 

(3)  The  ateady  etate  eolution  for  the  pitching  and  yawing  motion  la 
glvan  by  tna  third  vaotor  which  ia  rotating  at  the  rolling 
veloolty  of  the  mlealle.  Aa  a reault,  the  ateady  atata  roll 
orientation  of  the  mlealle  *o  the  lnatentaneoua  plane  of  pitch 
and  yaw  ia  a eonatant* 

(6)  For  atatlcally  atabl#  mlaailea  the  free  flight  motlonn  reaonate 
when  the  rolling  veloolty  of  the  mlaaile  end  the  nutation  rate 
approach  colnoidenoe. 

(7)  The  else  of  the  motion*  at  rtaontnoa  art  limited  only  by  tha 
dagret  of  dynaalo  atabl 11 ty  of  tha  mlaaile. 

The  ebove  characteristic*  of  the  eolutiona  for  the  frta  flight  motion* 


appear  to  form  a haul*  for  e eatlafeotory  explanation  for  the  three  phenomena 
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mentioned  la  the  INTROOUCTICN.  The  effect  of  the  rolling  velocity  on  the 
dispersion  of  the  missile  is  apparent  from  consideration  of  Eq.  (77)  snd 
Eq.  (73).  When  the  rolling  velocity  of  the  missile  is  zero,  large  dispersion 
may  result  from  the  quadratic  term;  however,  when  roll  is  introduced,  the 
quadratic  term  disappears  and  the  dispersion  may  be  reduced,  provided,  of 
course,  that  the  resonance  region  is  avoided,  Xt  is  this  resonance  condition 
that  appears  to  provide  an  explanation  for  the  second  and  third  phenomena 
mentioned,  since  the  large  angles  of  attack  and  yaw  and  the  resultant  trans- 
verse displacement,  predicted  by  the  theory  when  the  roll  rate  and  nutation 
rate  are  coincident,  tend  to  bear  out  the  observed  phenomena. 

Although  the  Tricyclic  theory  appears  to  contain  the  seeds  for  a satis- 
factory explanation  of  the  observed  physical  phenomena,  no  general  acceptance 
may  be  expected  until  the  ability  of  the  theory  to  represent  accurate  experi- 
mental data  of  free  flight  motion  has  been  thoroughly  investigated. 

Accordingly,  a program  haa  been  initiated  in  the  Spark  Photography1?” 20 > 3 3 
Flanges  to  obtain  the  required  accurate  free  flight  data.  The  general  program 
consists  of  an  investigation  of  motions  over  a large  range  of  rolling  veloci- 
ties and  includes  a detailed  investigation  of  the  resonance  region.  However, 
only  the  results  for  the  small  rolling  velocity  case  are  now  available  and 
are  given  in  the  following  section. 

EXPERIMENTAL  PROCSRAM 

The  purposes  of  the  present  preliminary  experimental  program  are  (l) 
to  investigate  the  ability  of  the  tricylic  theory  to  represent  the  actual 
free  flight  motion  of  missiles  having  Blight  configurational  asymmetry,  and 


(2)  to  determine  the  static  and  dynamic  aerodynamic  foree  and  moment  derive* 
tives  vhlcb  are  aeeociated  with  the  motion.  In  order  to  obtain  the  required 
free  flight  data  two  models  were  tested  in  the  Aberdeen  Spark  Photography 
Range.  17-80  (See  Pig.  13) 

The  models  employed  In  the  teats  had  a simple  arrov  configuration  with 
a cone-cylinder  body  of  fineness  ratio  of  10  and  with  cruel formed  single  wedge 
rectangular  fins  of  aspect  ratio  3 and  thickness.  One  set  of  fins  had  an 
angle  of  incidence  of  .3°  and  the  other  set  of  fins  had  no  incidence. 

(See  Pig.  l4).  The  models  were  launched  from  a special  railed  gun  which 
enables  launching  of  winged  and/or  finned  missiles.  (See  Fig.  13) 

RESULTS 

Motion 

The  experimental  data  for  the  pitching  and  yawing  motion  and  also  for 
the  transverse  displacement  of  the  models  as  obtained  from  the  shadowgraphs 
(Fig.  1 6)  taken  In  the  Range  are  given  by  the  points  in  Figs.  17,  lfl,  19, 
and  20.  In  these  same  figures  are  plotted  the  theoretical  curves  of  the 
tricyclic  motion  (Eqa,  (40)  and  (77)  which  were  "fitted"  to  the  experimental 
data  by  the  reduction  technique.**  30,32  ^ "fit"  of  the  theory  to  the 
experimental  data*  Is  given  by  the  Probable  Error  in  Table  II 
* The  Method  of  Differential  Corrections  is  used  to  "fit"  Eqs.  (40)  and  (77)  to 
the  experimental  data.  The  author  is  indebted  to  Mr.  C.  H.  Murphy,  Mathe- 
matician, B.R.L.,  for  programming  the  reduction  technique  for  automatic  computa- 
tion by  the  Bell  Relay  Computer,  and  also  for  suggesting  the  inclusion  of  the 
^ termB  in  EqB.  (21)  and  (27). 

**  See  Appendix  C. 
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FIG.  17-  EXPERIMENTAL  AND  THEORETICAL  FREE 
FLIGHT  PITCHING  AND  YAWING  MOTION, 
(RD.  2950) 
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♦ 01  radians 


FIG.  18  - EXPERIMENTAL  AND  THEORETICAL 

FREE  FLIGHT  PITCHING  AND  YAWING 
MOTION.  (RD.  2951) 


y 

INCHES 

FIG.  £0- EXPERIMENTAL  AND  THEORETICAL 
FREE  FLIGHT  TRANSVERSE  DIS- 
PLACEMENT. (RD.  2951) 


TABLE  II 


Thus,  since  the  Probable  Errors  of  the  residuals  is  of  the  same  order  of 
magnitude  as  the  estimated  error  in  measurement  in  the  Range  (Fig.  13)  the 
trioyclia  theory  of  pitching  and  yawing  motion  and  transverse  displacement 
may  be  considered  to  accurately  represent  the  actual  free  flight  pitching 
and  yawing  motions  of  these  models. 

Aerodynamic  Derivatives 

The  values  for  the  aerodynamic  derivatives  as  obtained  from  the  conetentB 
of  the  motion  are  given  in  Table  III  together  with  their  probable  error  aa 
obtained  from  the  least  squares  fit  of  the  theory  to  the  experimental  data. 


TABLE  III 


It  is  noted  that  the  values  for  the  moment  derivatives  due  to  angle  of 
attack  and  yaw,  C , as  obtained  from  the  motions  of  each  of  the  Identical 
missiles  differ  by  that  the  values  for  the  sum  ■>?  the  moment  derivative 

due  to  pitching  and  yaving  velocity  and  the  moment  derivative  due  to  rate  of 
change  of  the  angle  of  attack  and  yav,  , for  the  missiles  differ 


by  2.0^,  and  that  the  values  for  the  moment  derivative  due  to  control  surface 
deflection  (aeymmetry),  C , differ  by  0.5$.  Although  the  sample  is  small 

MSt 

ths  results  suggest  a good  reproducibility  of  the  static  and  dynamic  aero- 
dynamic derivative*  aa  obtained  from  the  Spark  Photography  Range  technique. 
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ipparoix  a 

Derivation  of  Basic  Differential  Equations  of  Motion 


The  derivation  of  the  basic  differential  equations  of  motion  by  employing 
the  total  kinetic  energy  of  the  missile  (Eq.  6)  and  the  Lagrange  equation 
(Eq.  7)  on  eaoh  of  the  coordinates  of  the  dynamical  system  (P#  0>r» 
z)  is  given  below: 

(1)  For  e : The  Lagrange  equation  for  6 is  given  by 


at  ) je  ’ 

Substituting  the  expression  for  the  total  kinetie  energy  of  the 
missile  (Eq.  6)  and  performing  the  indicated  operations  yields. 


(11) 


^rl1  ® ] “ [-  Xx  8 * zx  V1  e 008  6 “ 1 008  ® *Ln  3 " Q© 

Now  assuming  that  8,  f » 8,  and  ^ are  mall  and  that  their 
squares  and  products  may  be  negleoted  yields, 


I 8 + Ix  i - Qe 


(3A) 


(2)  For  ¥ : The  Lagrange  equation  for  [ft  is  given  by 

± fU)  - • v 

dt  ^ a yk/  9 y r 

Substituting  Eq.  6 and  differentiating  yields 

["  JX  $ ^ 6 + lX  ^ 0 f + 1 C0*2  9 ^ " [0]  " Q f 

or 

-I  J cos  0 8 - Xx  sin  0 0*  + Ix  sjn2  6^  + Ix  y/2  sin  8 cos  8 6 

n •!  • • 

+ I cos  9 y/  - I ^ 2 cos  8 sin  6 8 - Q y, 


(UA) 


(5A) 


(6A) 


(li)  For  x t The  Lagrange  equation  for  x ia  given  by 


d 


3T 

w 


9T 

3F 


Substituting  Eq.  6 and  differentiating  yields 


& C-iJ  -Co] 


m x ■ Qx 

(5;  For  y r The  Lagrange  EquaV  ->n  for  y is  given  by 


Substituting  Eq,  6 and  differentiating  yields 

»y  - 
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The  Lagrange  equation  for  * ie  given  by 


3 T 
Ts 


- Q 


a 


Substituting  Sq.  6 and  differentiating  yields 


(l6i) 


- [o]  ■ 9, 

••  . 

» * - Q, 

Now  multiplying  Bq«  (7A)  by  i and  adding  to  Bq.  (3A)  yields 

i ( e ♦ i f ) - i ix  flf  ( e + 1 yo  • Qe  * 1 * P 


(17A) 

(18A) 


(19A) 


Introducing  oomplex  -rariables  by  defining 

Si  - i ♦ i f 
n - e + i 


and  substituting  into  Eg.  (1?A)  yields 

i/i  * 1 ix.  i si  - qa 

Also  multiplying  Eq.  (l8A)  by  i and  adding  to  Eq.  (l?A)  yields 
m (y  + i *)  " ^ + 1 


defining 

S ■ y + i * 

S ■ y + i » 

••  ••  • • 

S » y + i* 


(20A) 

(21A) 


(22A) 


(23A) 

<2UA> 

(2*A) 

(26A) 


:*yr.  tj  «ra  . 
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APPENDIX  B 

Aerodynamic  and  Ballistic  Nomenclature 

The  aerodynamic  force a and  moments  in  both  the  Aerodynamic  nomen- 
clature^ ami  the  Ballistic  nomenclature,  and  also  in  the  Aero-Ballistic 
nomenclature  (for  missiles  possessing  trigonal  or  greater  rotational 
symmetry  and  mirror  symmetry)  ore  given  in  Table  I. 

The  relation  between  the  Ballistic  J's  and  the  Aero-Ballistic  C's 
are  given  in  Table  II. 

From  the  definitions  and  relations  of  Tables*  I and  II,  some 
of  the  basic  equations  of  the  report  may  bo  converted  to  Ballistic  form. 

The  oonstants  of  the  Aerodynamic  force,  Eq.  (21),  and  the  aero- 
dynamio  moment,  Eq.  (27)*  in  Ballistic  nomenclature  are  given  by: 

FORCE  CONSTANTS 


a - - p u2  d2  + i Kr  p u w d3 

b - /f  XF  P » d*4  + i K a pud3 

^ LN  P 11  ^ d P W d 

d - • K LJF  p ^ " i/^LS  P ^ 

e • p u2  d2 

l 


rs 

s ' 


|’*v,v  s>Hopgood,  R,  C.,  "A  Proposed  Revision  of  American  Standard  Letter 

p-.-v  Symbols  for  Aeronautical  Sciences",  Aeronautical  Engineering  Review. 

anuary,  19?. 3. 
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TABLE  I 

AERODYNAMIC  FORCES 


aerodynamic 


Yg  P - Cy  r?s 

F p 

2c  “ ■ °2 

Yp  a p0  ■ cy  “qS 

pa 

ZpP  pP  ' 

Yr  r “ q s 

2,  1 ■ \ 

*pq  P<*  " CYb  $>$>  q S 
pq 

Vpr  ■ 1 

H p “ CY'  $>  q s 

zi  • ■ Oj, 
e 

\ pa  ■ cyp.  ®)(^)  q s 

Zpj§  Pf5  - 0 

Y-  * " °i£  5 s 

Zq  * " CZ. 

Q 

Ypq  P*  " Crp.  q s 

v ^ • °i 

\ Vcy  ^qs 

Tl 

\ h ' °j 

pP 


2bwrb, 


- 


Sp 


/Pb\  - 


pr 
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TAiiLti  I 


AERODYNAMIC  F CROSS 


AERO  - BALLISTIC 

BALLISTICS 

N$£  “ * 8 

N - - Kjjpu2  d2£ 

NPJ  pJ  “ 1 °Npfl  ’ S 

F**iKj,pucod^^ 

W 1 °Nq  <ty  < S 

S - i Ks  p u d3Y^ 

Nprf  CN_  ^H^)qS 

Pq 

XF  - Kjjj.  p o>  d1*  ^ 

m p u d3 1> 

y i ■ 1 %.  ##  5 5 

LF«-iKLFptodilg 

“+|A- 1 $)(h>  i ■ 

Q 

ls  - - i p dunri 

"p+fA"  °N  . * 8 

v pq 

- - Klxf  p S dH 

V«  ‘ \ ®s  5 s 

Ne  - \ p ”2  d2  6t 

TABLE  I 


«l 


AERODYNAMIC  MOMENTS 


c * Q.  a q S c 
a nc 

Np  P - °np  P q S t 

>W  pP  ■ \p(lr>  P 5 s « 

Npa  pa  " 8 « 8 b 

M,  « ■ \ # 5 s « 

Nr  * " %#  * 8 b 

v pr  ■ « s ° 

NPq  Pq  " Cnpq(57)(l5>  < 8 b 

»; 4 ■ cni  # 5 3 • 

u 

Np  f " °n^  * 8 b 

mpP  p*  “ ^ q « o 

NPi  P8-°V#^^Sb 

Mi  * " Cm-  (^)(W)  * S 0 

Nr  * " °n*  (?T)(W)  « 3 b 

V p*  " °v  5 s c 

Npq  pq  " CHp^  « 8 

\ 6E  ’ c»^  6B  * S 0 

\ " °n^  «fe  » 8 b 

WHERE  ^ ■ p + i a,  ^ ■ p + i a , 
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TABLE  I 


AERODYNAMIC  MOMENTS 


AERO-BALLISTIC 

BALLISTIC 

1 °Ma  ^ ^ S 0 

M - -i  Kj^  p u2  d3  Jl 

v*  ■ 5 5 0 

T - -Kj  p u u d14  £ 

V ■ <>Mq$>  ’ 8 0 

H«*KHpud^ 

Mpjj  P»  ■ -1^  $>$>  5<« 

XT  - -iKXT  p (o  d*>t 

»ii  ■ -ioh.  w « 8 « 

m - iKm  p u du  j 

M p£  ■ v $)(&  5 8 ° 

LT  ■ K^T  p « d*£ 

• °m.  « 8 ° ' 

<1 

LH  ■ Kjjj  p d*  V( 

% " -1CMp4  (!*)<R)<W>  « 8 0 

LXT  - -Ujjj  Pgd6^ 

V eA  - i V fl4  q fl  o 

Me  - % p *2  d3  6C 

ILI 


MOMENT  CONSTANTS 

k ' 2 ^ 

A ■ - KT  pu^d4  • i ^ P«  d 

ll  *9 

B - - Kh  pud4  + i KjT  p A)  o 

C - Kjj,  p«dS  ♦ iK^p  udlu  (2B) 

t 6 

D - ^ P 4s  • 1 Im  P «'u  1 


- i 


H 


p»2  d3 


Substituting  tbeae  expressions  to  Eqs.  (32)»  «nd  (3U)»  and  eorabining 
yields  the  differential  equation  far  the  pitching  and  yawing  motion 


• 0 

$ 


V 

- T 


j - - k"2  ^JH  ♦ Jjjj)  + 1 V D'2  (Jlt  + JxT)  + JF  * A^j 

f-k”2 


-2 


V 


A/b  Jf  * i > A/p  J«  ♦ 1 v k Ji 


APPENDIX  C 


Reduction  Technique 

It  is  the  purpose  of  the  reduction  technique  to  "fit"  the  Tricyclic 
Theory  to  the  experimental  data  obtained  from  the  Aerodynamics  Range  and 
to  determine  the  various  constants  associated  with  the  motion.  The 
present  procedure  is  to  first  consider  the  pitching  and  yawing  motion 
and  than,  using  the  values  of  the  constants  obtained  from  the  pitching 
and  yawing  motion,  reduce  the  transverse  displacement.  Once  these  con- 
stants have  been  obtained  the  representation  of  the  experimental  data  by 
the  theoxy  (i.e.,  "fit")  way  be  determined  and  the  aerodynamic  deriva- 
tives may  be  calculated. 


Yawing  Motion 

For  the  purpose  of  Range  reduction  the  equation  of  yawing  motion 
(Rq.  i|0)  is  modified  in  two  weysi 

(1)  the  independent  variable  is  changed  from  time  to  distance 
along  the  trajectory  of  the  missile,  a. 

(2)  the  pitch  and  yaw  are  measured  from  the  instantaneous  velocity 
vector  of  the  center  of  gravity  to  the  miseile  and  are  taken 
as  positive  in  *he  first  quadrant  looking  at  the  approaching 
missile*,  and  is  written  ae 

„ „ *•**(•*,  * V *>  ifc  ♦ b,  .)  ■)  ik  ♦ b. 

p+io-e  1 1 e ^ 1 ♦ e 2 e 2 2 


2. 


+ e 


3 * V 


) 


(1C) 


where  the  new  constants  are  related  to  those  in  Bq.  (U6)  by 


io«io  hi 

^l 

\ " r.w? 

lo®10  1 ^2  1 

*•2 

\ " r55T? 

lo«io  M 

bi  " r (2°H5c)(7C) 
b2  - Jr  (30) (60) (80) 
b3  ■ $ (U0)(9C) 


* In  the  standard  system  the  yaw  is  measured  from  the  missile  to  the 
instantaneous  velocity  vector  of  the  center  of  gravity  and  is  taken 

as  positive  in  the  third  quadrant  looking  at  the  rear  of  the  missile 
(See  Fig.  3). 
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Substituting  the  Range  data  for  p,  a and  z at  each  station  into  Eq.  (1C) 
yields  a set  of  5 0 equations  which  are  non-linear  in  the  constants  to  be 
determined  thus  the  Method  of  Least  Squares 3^  cannot  be  applied. 


Differential  Corrections 


31 


Eq.  (1C)  is  therefore  expanded  in  a Taylor's  Series  in  which  the 
higher  order  terms  are  heglacted  as 

2-302(*k1  * \ ‘)  iL  ♦ b,  «). 

‘3m -lo-*  ° ° • ® ° [2.302^  ♦ 2.302^ 

+ + i zAb^  | 

2,302(\  * v “)  2l\  * \ >)r 

* e 0 e 0 0 ]2.yi2A^  + 2,302^^ 

+ idag  + i z^b2  1 

2*302M  i/e«  ♦ b«  s\ 

f a 0 e 0 0 J2.302  *K  * 1 *3  ♦ i * b^J 

where  £ ^ ■ measured  value  of  angle  attack  and  yaw  in  Range 

% 0 ■ value  computed  from  Eq.  (1C)  using  initial  values  of 
constants 


(10c) 


Since  these  equations  *re  linear  in  the  differential  corrections  the  Method 
of  Least  Scares  may  be  applied  for  their  determination,  provided  that 
estimates  0 i initial  values  of  the  constants  may  be  made.  Once  valuee  of 
the  differential  corrections  are  determined,  they  may  be  added  to  the 
original  initial  values  of  the  constants  and  this  process  of  Differential 
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it  1b  noted  that  initial  estimates  of  thas®  constants  nay  be  readily 
determined  for 

(1)  "by  the  rolling  velocity  of  the  miaaile  as  a function  of  t 

may  be  measured  directly  from  the  observed  motion  in  the 
Range1?*20 

(2)  a.y  the  angular  position  of  at  z may  be  obtained  from  a 

knowledge  of  t he  roll  orientation  of  the  missile  at  a ( see 
above  (1)  and  a knowledge  of  the  orientation  of  K,  with 

respect  to  the  asymmetry*  vhioh  is  fixed  in  the  missile  (kkj.T!?) 

(3)  e2*302  l\)  i the  size  of  at  zQ  may  be  estimated  from  Eq.  (U3) 

provided  that  the  asymmetry  is  known  from  the  physical  measurements 
of  the  missile  and  th4  estimated  values  of  toe  aerodynamic 
coefficients  are  available.  In  maty  oases,  inspection  of  the 
experimental  pitching  and  yawing  motion  yields  a good  indication  of 
the  size  of  Kj. 

Thus  initial  values  far  tu  , a,  and  b,  may  be  obtained.  Subtracting  the 

»» .5  22 

third  term  from  the  experimental  data  and  applying  the  standard  technique''* ■ “ 
yields  estimates  of  the  remaining  initial  values  of  the  constants. 

Transverse  Displacement 

The  reduction  of  the  transverse  displacement  of  the  missile  follows 
directly  from  tho  reduction  of  the  pitching  and  yawing  motion.  Writing  I2q. 
(77)  as 


(Vv  2 (jWz  (sVv)z 

S *>  Kj  e * k2  e ♦ kj  e ♦(kj|/v)jT  ♦ kj  13C) 


and  substituting  the  .lan-e  data  for  a and  Z end  the  pitch  and  yaw  contents 
previously  doterrv*  led  yields  a set  of  Z ? equations  which  are  lino?.:’ 
in  tho  unknowns,  k,,  k?,  k_,  ki  and  k^jano  thus  the  Method  of  Least  squares 
may  be  applied  f.-.-Hheir  determinations. 


"The  angular  orientation  of  the  asymmetry  with  respect  to  the  missile  may 
be  determined  by  measurement  prior  to  firing. 
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Corrections  successively  repeated  until  the  sum  of  the  squares  of  the 
residuals  is  a minimum.  The  values  of  the  final  constants  so  obtained 
may  be  considered  as  the  best  the  data  can  yield  employing  the  theory 
of  motion.  If  the  probable  Error  from  the  residuals, 


PE  *, 


(11C) 


is  of  the  same  order  of  magnitude  as  the  estimated  error  in  measurement 
of  the  experimental  data,  then  the  theory  may  be  considered  to  represent 
accurately  the  data. 


Initial  Values 


Hie  determination  of  the  original  initial  values  of  the  constants 
is  critical  for  employing  the  Method  of  Differential  Corrections^  for 
if  their  determination  is  poor  then  the  process  may  not  be  convergent.* 

The  technique  used  in  obtaining  these  initial  values  is  to  remove 
the  effect  of  the  third  term  from  the  experimental  dataj  then  the 
resultant  data  will  be  epicyclic  and  the  standai  l technique  for  deter- 
mining the  initial  values  may  be  used.^S2 

The  effect  of  the  third  arm  on  the  experimental  data  may  be  sub- 
tracted once  the  constants  of  this  term  are  estimated. 

Considering  this  tem 

2.3°2(.  ) i(i  . b ,)  j-  -I 

« J.  ° 0 1^.3024.*  .n.j  « 1 .AbjJ  (120) 

where  e * ■ size  of  arm  Kj 


• angular  position  of  arm  at  z0 

b,  - rolling  velocity  of  the  missile  as  a function 
J of  i 


* Of  course  if  the  theory  is  not  correot,  the  Process  may  also  be  di- 
vergent. 
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